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THE OXYGEN ABUNDANCE IN THE INNER H II REGIONS OF M101. 
IMPLICATIONS FOR THE CALIBRATION OF STRONG-LINE METALLICITY INDICATORS 1 

Fabio Bresolin 

Institute for Astronomy, 2680 Woodlawn Drive, Honolulu, HI 96822; bresolin@ifa.hawaii.edu 

ABSTRACT 

I present deep spectroscopy of four Hll regions in the inner, metal-rich zone of the spiral galaxy M101 
obtained with the LRIS spectrograph at the Keck telescope. From the analysis of the collisionally excited lines 
in two of the target Hll regions, H1013 and H493, I have obtained oxygen abundances 12 + log(0/H) = 8.52 
and 12 + log(0/H) = 8.74, respectively. These measurements extend the determination of the oxygen abundance 
gradient of M101 via the direct method to only 3 kpc from the center. The intensity of the C II A4267 line in 
H1013 leads to a carbon abundance 12 + log(C/H) = 8.66, corresponding to nearly twice the solar value. From 
a comparison of the continuum temperature derived from the Balmer discontinuity, T(Bac) = 5000 K, and the 
line temperature derived from [O III] A4363/ A5007, T[0 III] = 7700 K, an average temperature T = 5500 K and 
a mean square temperature fluctuation t 2 = 0.06 have been derived. Accounting for the spatial inhomogeneity 
in temperature raises the oxygen abundance obtained from the oxygen auroral lines to 12 + log(0/H) = 8.93. 
These findings are discussed in the context of the calibration of strong-line metallicity indicators, in particular 
of the upper branch of /?23- There is no evidence for the strong abundance biases arising from temperature 
gradients predicted theoretically for metal-rich H II regions. 

Subject headings: ISM: abundances — galaxies: ISM — galaxies: abundances — galaxies: individual (M101) 



1. INTRODUCTION 

As probes of the gas-phase chemical composition of star- 
forming galaxies, Hll regions are crucial targets in the 
quest to understand the chemical evolution of the Universe. 
Through the analysis of bright H II regions in external spiral 
and irregular galaxies it is, in principle, rather straightforward 
to measure the abundances of elements such as oxygen, nitro- 
gen, neon, argon and sulphur, and therefore constrain models 
of nucleosynthesis in massive stars and of galactic chemical 
evolution. 

In recent years a great deal of interest has been raised con- 
cerning the derivation of chemical abundances in metal-rich 
H II regions, i.e. those where the oxygen abundance reaches 
or exceeds that measured in the Sun [12 + log(O/H) = 
8.66, lAsplund etal] 1200411 2 . This work is motivated by the 
importance of obtaining accurate abundances for the mea- 
surement and the interpretatio n of galactic chemical abun- 
dance gradients (Pilvugi n et alJl2004l |Chiappini et al .1 120031 
Carigi et al.l 120051) . and for constraining the metal-rich end 
of the luminosity-metallicity (| Kobulnickv & Kewlev |2004 t 
ISalzer et alJl2005t iMaier et al.ll2005l iLamareille et al.ll2006l) 
and mass-metallicity j Tremonti et al. 2004; Savagli o et al.l 
2005; lErb et all 120061) relations of galaxies up to red- 
shifts z ~ 2. A large, solar-like oxygen content is in- 
ferred for star-forming ga laxies both locally (BresolinetaL 
| 2004t fPilvugin et al. 2006) and at high redshift ( Shanle v et al. 
2004). In addition, e mpirical metallicity trends in the Wolf- 
Rayet WC/WN ratio (Mevnet & Maeder 2005; Masse\ l2003l) 
and in the Cepheid Period-Luminosity Relation dSakai et aD 
2004) depend critically on the reliability of the abundances of 

1 The data presented herein were obtained at the W.M. Keck Observatory, 
which is operated as a scientific partnership among the California Institute of 
Technology, the University of California and the National Aeronautics and 
Space Administration. The Observatory was made possible by the generous 
financial support of the W.M. Keck Foundation. 

The terms metallicity and oxygen abundance will be used interchange- 
ably throughout the paper. Abundances are measured relative to hydrogen. 



metal-rich H II regions. 

The difficulty of measuring chemical abundances from 
the classic analysis of optical forbidden lines (collision- 
ally excited) in high-metallicity H II regions derives from 
the increased proportion of gas coolants, in particular oxy- 
gen, which emits line radiation mostly via the fine-structure 
[Olll] AA52, 88 fim lines. As a consequence of the re- 
duced gas temperature, the auroral-to-nebular line ratios 
[e.g. [O III] A4363/( A4959 + A5007)], employed to determine 
the electron temperature T e , become progressively smaller 
with increasing metallicity, due to the exponential dependence 
of the collisionally excited line emissivity on T e . Gener- 
ally, at metallicities approaching solar the [O III] A4363 au- 
roral line is too faint to be observed in extragalactic nebulae 
even with lOm-class telescopes, as a result of the strong T e 
dependence, and of the fact that the far-IR [O III] lines be- 
come stronger at the expense of the optical lines. However, 
auroral lines from other ionic species can be detected in the 
optical spectra of H II regions at relatively high abundances, 
such as [S III] A63 12 and [Nll]A5755, as a consequence of 
their larger Boltzmann factors exp(~x/kT e ) (x is the energy 
difference, a few eV, between the two levels involved in an 
auroral line transition). Still, the detection of these diagnostic 
lines in extragalactic nebulae requires high surface brightness 
and, generally, lOm-class telescope capabilities. 

Wit h the ai d of photoionization models (Stasihska 
119781 iGarnettl 119921) . one can relate the electron 
temperature derived for the [S III] or [N II] emit- 
ting regions [from [S III] A6312/( A9069+ A9532) and 
[Nll]A5755/(A6548+ A6583), respectively] to the tem- 
perature of the [O III] emitting region. When A4363 is 
unavailable, this allows the measurement of the oxygen 
abundance from the intensity of the strongest lines only 
(AA4959, 5007), having fixed the line emissivity from the 
knowledge of the electro n temperature . This techn ique has 
been recently applied by Bresolin et al. (2004, 2005) in Keck 
and VLT studies of the chemical abundances of samples of 
extragalactic metal-rich H II regions. 
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The availability of direct abundances is essential for 
the empirical calibration of the so-called strong-line 
methods, in which the strength of easily observed neb- 
ular emission lines are related to the metallicity of the 
emitting regions. Among the various methods proposed, 
I recall jfo = (TO ill A37 27 + [O m] AA4959, 5007)/H/3 
jPagel et alJ fl979h . N2 = [ Nlll A6583/Ha 

JstorchTBergmann et alJ Il994t iDenicolo et alJ 

120021). 03 N2 = ([ in] A5007/H/3)/(rNi il A6583/Ha) 

jAlloin et alJ fl979[ iPettini & Pagel l200l . 

59C( = (rSHlAA6717 673 1 + [S in] AA906 9, 9532)/H/3 

(Vflchez & Estebanl 119961 iDiaz & Perez-Monteroj 
2000), the P-method jPilvug id l2000t iPilvugin & Thuanl 
120051) and Ar 3 3 = [A r ml A7135/rOiii l A5007, 

S3O3 = [S in] A9069/[O in] A5007 ( Stasihska! 120061) . For 
in-depth discussions of these abundance indi cators the reader 
is referr ed to these papers, as w ell as to IKewlev &D opita 
(2002), iPerez-Montero & Dfazl J2005I) and iDopita et al] 
(2006a). The strong-line methods allow the estimation 
of nebular chemical abundances in cases where a direct 
approach, based on the detection of the auroral lines, is not 
feasible, and in general when the signal-to-noise ratio of the 
spectra is sufficient to measure only a restricted number of 
emission lines. This approach is widely adopte d for the deter- 
minat i on of abundance gradients in galaxies ( P ilvugin et alJ 
120041 120061) and of oxygen abundan ces in in termediate- 
and high-redshi ft star-forming galaxi es ( Shaplev et al] l2004t 
ISavaglio et aDl2005t iMouhcine et alJ 120061) . The popularity 
of these metallicity indicators is also due to the increasing 
necessity of estimating chemical abundances for galaxies 
observed as part of large s pectroscopic surveys, such as the 
Sloan Digital Sky Su rvey ilTremonti et alJl2004t iLiang et alJ 
2006; ShiSMi). 

Calibrations of strong-line metho ds based on grids 
of photoio nizati on models ("e.g. [Edmunds & Pagel l Il984t 
McGa ughUl991t IZaritskv et alJ 11994 ICharlot & Lon ghetti 
120011 IKew lev & Dopita 2002) provide abundances that 
are, at the metal-rich end, considerably larger (up to 
0.5 dex) than those obtained from direct measurements 
JCastellanos et al J 120021 IKennicutt et al]l2003l IGarnett et alJ 
2004b: iBresolin et alJl2004t 120051) or from ^-calibrated em- 
pirical methods. Mixing T^-based abundances (routinely done 
at low metallicity) with abundances from model grids (at 
large metallicity) leads then to spurious discontinuities in 
abundance diagnostic diagrams or breaks in the luminosity - 
metallicity relation (Nagao et al. 2006; Lee et al. 2006J). The 
reasons for the discrepancy are still currently poorly known: 
additional sources of heating need to be postulated in or- 
der to re produce theoretically the measured \Q ill] A4363 li ne 
strength JStasinska & Schaereil199ftlT~uridiana et al. 1999). 

There are a number of open issues regarding the study of 
metal-rich H II regions that motivate further investigation of 
ionized nebulae in nearby galaxies. First of all is the assess- 
ment of the importance of the effects of temperature gradients 
in H II regions and the resulting abundance biases, expected 
from theoretical models (Stasihska 1980). As the metallicity 
approaches the solar value, strong temperature gradients de- 
velop in H II regions, due to the more efficient cooling in the 
central zone from far-IR [O III] fine-structure lines relative to 
the outer region. Since the line strength is weigthed towards 
warmer regions, due to the strong T e dependence of the line 
emissivities, the net effect of these temperature gradients is 
that the electron temperature derived from the auroral lines 
over-estimates the real temperature, and consequently the 



oxygen abundance is systematically under-estimated (IGarnett! 
1992; Stasihska 2005). In general, in the presence of spatial 
temperature variations within H II regions, methods based on 
collisionally excited lines can lead to under- estimate chemi- 
cal abundances by a factor of two or more dPeimbertlll967t 
iTorres-Peimbert et alll980l) . 

Alternative methods of abundance determination, less de- 
pendent on the nebular temperature structure than collision- 
ally excited lines, need therefore to be explored in extra- 
galactic H II regions of high metal content. Metal recombi- 
nation lines and infrared fine-structure lines offer such alter- 
natives. In both cases the line emissivity is only moderately 
dependent on T e . There are obvious difficulties in exploiting 
these techniques, however: recombination lines are very weak 
and h ard to detect in extragalactic H II regions (Esteban et al. 
2002) , and infrared observations require orbiting telescopes 
JGarnett et alJ2004al) . 

Studies of the Oil recombination lines in Galac- 
tic and extragalactic Hll regions find oxygen abun- 
dances 2 to 3 times larger than those obtained from 
[Olll]A4363/(A4959+ A5007), and conclude that the dis- 
crepancy can be e xplained by the pre s ence of temper- 
ature fluctuations ( Peimbert et al] 119931 iPeimbertl 120031: 
iGarcfa-Roias et alJ2006l) . There is still no general consensus 
about the source of temperature fluctuations in Hl l regions 
and planetary nebulae (Peimbert & Peimbert] 120061) . and al- 
ternative explanations for the abundance discrepancy between 
collisionally excited lines and recom bination lines ha ve been 
proposed tLiu et al.l2000tlTsamis & Pequignol2b 05). 

In this paper I report on new deep spectroscopic obser- 
vations of four Hll regions in the inner zone of the spi- 
ral galaxy M101. The radial oxygen abundance gradient in 
this galaxy has been studied in detail by previou s authors 
iTorres-Peimbert et alJl!989t IKennicutt et alJl2003i hereafter 
KBG03). The oxygen abundance of the innermost H II re- 
gions studied thus far, using auroral line detections, is ap- 
proximately 12 + log(0/H) = 8.7, with an extrapolated central 
abundance 12 + log(0/H) = 8.76 (KBG03), corresponding to 
1.26(O/H) . The new observations presented in this paper 
provide additional high-quality spectroscopic data for some 
of the most metal-rich Hll regions in this galaxy, with the 
goals of: (a) obtaining direct abundances closer to the galaxy 
nucleus from auroral lines; (b) verifying, if feasible, these 
abundances using metal recombination lines, and (c) mea- 
sure the effects of temperature fluctuations on the direct abun- 
dances. Observations and data reduction are explained in §2. 
The abundances derived from collisionally excited lines in 
two H II regions are presented in §3, while those derived from 
recombination lines for helium and carbon for one of them 
are presented in §4 and §5, respectively. A measure of the 
temperature fluctuations in one of these H II regions, and their 
effect on the chemical abundances derived from collisionally 
excited lines, is given in §6. In §7 1 discuss the implications of 
the results obtained for the calibration of strong-line metallic- 
ity indicators. The main results of this paper are summarized 
in §8. 

2. OBSERVATIONS 

The four targets chosen for this investigation are among the 
brightest H II regions in th e central 3 arcm inutes (~ 6 kpc 
at a distance of 6.85 Mpc, Freed man et alJl2001l) of M101: 
H493, H507 (also Searle 2. from ISearlelll971l) . H972 and 
H1013 (Searle 3). Tabled contains the identificat ion and the 
celestial coordinates taken from the catalog of Hod ge et"aT1 
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( 1990), together with the deprojected galactocentric distance 
in units of the disk i sophotal radius (Rq = 14'A), taken from 
iKennicutt & Garnettl (119961) . 

The spectra of the Hll regions were obtained on May 12, 
2005 with the Keck I telescope on Mauna Kea using the Low 
Resolution Imaging Spectrometer (Oke et al. 119951) . The ob- 
serving conditions were photometric, and the seeing was sta- 
ble at O'.'l throughout the night. Using a dichroic splitter, blue 
and red spectra were obtained simultaneously over most of the 
optical wavelength range. For the blue spectra (3300-5600 A) 
a 600 lines mm -1 grism blazed at 4000 A was employed, 
yielding a spectral resolution of ~5 A with a 1.5 arcsec-wide 
slit. In the red (4960-6670 A) a 900 lines mm" 1 grating blazed 
at 5500 A was used (4 A resolution). In order to extend 
the spectral coverage to the near-infrared [S III] lines, spec- 
tra were obtained also with a 400 lines mm -1 grating blazed 
at 8900 A, covering the 6050-9800 A wavelength range at 9 A 
resolution. 

The total integration time for the blue and red spectra 
ranged between 4200 and 4500 seconds, split into two or three 
exposures. For the near-IR spectra single 300 s or 900 s expo- 
sures were obtained. I also acquired 300 s blue and red spec- 
tra for the two brightest nebulae, H1013 and H972, in order to 
avoid saturation of the Ha emission line. Table^summarizes 
the exposures secured for this program. 

The airmass of M101 varied during the observations be- 
tween 1.21 and 1.83. In order to minimize the effects of dif- 
ferential atmospheric refraction the position angle of the slit 
was adjusted to match the parallactic angle as it varied during 
the course of the night. Observations of the standard stars 
LTT 7987, BD+28d4211 and BD+25d4655 were obtained 
during twilight for flux calibration. 

2.1. Data reduction and line flux measurement 

The long-slit data reduction was carried out with IRAF 3 , 
version 2.12, using the PyRAF 4 command language. The 
spectra were corrected for bias, flat-fielded and wavelength 
calibrated using Hg + Zn + Cd lamps (blue range) and Ne + Ar 
lamps (red range ). For the flux calibratio n the Mauna Kea ex- 
tinction curve of Krisciunas et al. ( 1987) was employed. Due 
to the lack of proper calibration data, the flux calibration at 
wavelengths longer than 9200 A was not carried out. Spec- 
tral regions in common between adjacent wavelength ranges 
showed a good agreement in the calibrated continuum level 
and in the flux of lines in common, indicating a 5% accuracy 
in the relative flux calibration. 

The strength of the emission lines was measured by in- 
tegration of the flux under the line profile, measured be- 
tween two continuum points selected interactively. These 
fluxes were cor rected for interstellar reddening adopting the 
iHowarthl (119831) analytical formulation of the iSeatonl (119791) 
law. The reddening coefficient C(H/J) was obtained from the 
Balmer decrement, considering the measured intensities of 
Ha and H ^ relative to H/3, compar ed to the case B values 
taken from Storev & Hummer ( 1995) for an electron temper- 
ature T e = 7500 K. The equivalent width of the Balmer lines 
in absorption originating from the underlying stellar compo- 

3 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 

4 PyRAF is a product of the Space Telescope Science Institute, which is 
operated by AURA for NASA. 



nent was adjusted iteratively in order to reach agreement in the 
C(H/3) determined separately from the Ha/H/3 and H7/H/3 ra- 
tios. Higher order lines of the Balmer series are progressively 
more sensitive to the amount of underlying absorption, and 
were thus not included in this procedure. 

The result of the line flux measurements, corrected for red- 
dening, is presented in Table [2] The quoted uncertainties ac- 
count for statistical errors, as well as the uncertainty in the flux 
calibration, the flat-fielding, and the reddening. I have com- 
pared the flu xes of the strongest metal lines with the values 
published by Kennicutt & Garnett ( 1996). In general there is 
agreement at the 1-er level. The comparison is slightly worse 
(2-cr) for the [0 11] A 3727 line in H1013 a nd H507. In simi- 
lar comparisons (e.g. Bresolin et al. 2005) the [Oil] A3727 is 
often found to be the one showing the largest discrepancy. It 
is not clear if this is due simply to errors in the line flux mea- 
surements and/or reddening corrections, or if variations in the 
regions sampled by the various slits can also account for the 
observed differences. 

As Table|2]shows, auroral lines have been detected in three 
of the four targets: H493, H507 and H1013. These lines 
will be used in Section 3 to derive electron temperatures and 
ionic abundances from the collisionally excited lines. The 
[N II] A5755, [S III] A63 12 and [0 11] A7325 auroral lines have 
already been measured in H1013 by KBG03. For these lines 
the comparison with the fluxes presented here is excellent. 

The spectrum of H507 shows a large number of Fe II lines 
(Fig. [0, resembling the emission-line rich spectra of known 
luminous Galactic iron stars, such as HD 316285 and ri 

J^^^^^^^^^^^~ ^^^^^^^^^^^^ ^^T^ 

Carinae (Walborn & Fitzcatrick 2000). These are massive 
stars believed to be in an advanced stage of evolution. The 
presence of evolved massive stars in H507 is also evident 
from the Wolf-Rayet emission features at ^4650 A (both 
the N III AA4634-42 and the He II A4686 are detected, indica- 
tive of WN stars) and 5696 A from cm (late WC stars). 
In fact, W-R emission lines have been detected in all four 
Hll regions studied here. The high detection fraction of W- 
R features, often including the signatures of WC stars, in 
high-metallicity H II region s has been noted in earlier works 
(Bres olin et alJl2004l 120051 lHadfield et aUl2005l) . and is ex- 
pected frojnJhe_jncreasem_the stellar mass loss with metal- 
licity (Mevnet & Maeder 2005). Given the scope of this pa- 
per, I will not further discuss the W-R content of the H II re- 
gions under consideration, and I will simply note that the stel- 
lar envelope of stars with Fe II lines in emission can also be 
responsible for e mission in auroral lines such as [N II] A5755 
and [ O Hi A7325 llHilher et al.ll 19981 iBorges Fernandes et all 
120011) . Therefore, these lines will not be used to measure neb- 
ular electron temperatures in H507. 

Electron densities were estimated from the 
[S II] A6716/A6731 line ratio, which in all Hll regions 
is consistent with the low-density limit. This is also the 
case for the [CI m] A5518/A5538 ratio measur ed forH1013 
comp ared to the theoretical low-density value ( Keenan et alJ 
2000). A common density N e = 100 cm -2 will therefore be 
adopted in the remainder of the paper for all H II regions. 

3. METAL ABUNDANCES FROM COLLISIONALLY EXCITED LINES 

The following analysis focuses on the two Hll regions 
for which reliable electron temperatures could be measured: 
H493 and H1013. For the derivation of ionic abundances 
from c ollisionally excited lines I relied on the programs by 
IShaw & Dufourl(ll995l) available in the STSDAS nebular 'pack- 
age running under IRAF. The atomic parameters and the col- 
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TABLE 1 

Object sample and exposures acquired 





R.A. 


DEC. 






Exposures 




Object 


(J2000.0) 


(J2000.0) 


R/Ro 


3300-5600 A 


4960-6670 A 


6050-9800 A 


H493 


14:03:03.4 


+54:21:25 


0.10 


3 x 1500 s 


3 x 1500 s 


1x900 s 


H507 (Searle 2) 


14:03:04.2 


+54:19:28 


0.14 


3 x 1500 s 


3 x 1500 s 


1x900 s 


H972 


14:03:27.8 


+54:21:31 


0.16 


2x2100s 


2x2100s 


1x300 s 










1x300 s 


1x300s 




H1013 (Searle 3)... 


14:03:30.7 


+54:21:14 


0.19 


3 x 1500 s 


3 x 1500 s 


1x300s 










1x300 s 


1x300s 
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FIG. 1 . — The spectrum of H507 contains a large number of Fe II lines, some are shown here in the 4150-4830 A (top) and 5150-5830 A (bottom) wavelength 
ranges. The Wolf-Rayet features at ~ 4650 A and 5696 A are also labeled. 



lision strengths adopted are those implemented in the 1997 
version of the software, except for the S HI col lision strengths, 
which were taken from lTaval & Ouptal (fl999>> . 

In the case of H1013 five values of the electron tempera- 
ture were derived from the ratios of the available auroral lines 
to the corresponding nebular lines (see Table [3}. The values 
from the different ions are all similar to each other, and are in 
the 7680-8360 K range. The ionic abundances were derived 
adopting T[0 III] for ++ and Ne ++ , T[S III] for S ++ and Ar++, 
T[S II] for S + , and a weighted average of T[0 II] and T[N II], 
corresponding to 8070 ± 120 K, for + and N + . 

For H493 only T[S III] and J [Nil] could be directly mea- 
sured. In this case, therefore, I determined the temperature 
relative to the emission zone (assumed to be equal to the 
temperature measured from th e [O Hi] lines) followin g previ- 
ous papers in this series ( Bresolin et al. 2004|, 12005ft . i.e. by 
adopting the relationship between 770111] and TTSlll] pre- 
dicted by photoionization models dGarnettlll992l Istasiriska 



QUI: 



T[Sm] =0.83 T[Om] + 1700 K 



(1) 



where the electron temperatures are in K. However, because 
oxygen is present in H493 almost exclusively as + , this pro- 
cedure has a very little impact on the derived total oxygen 
abundance. The other assumption made in the abundance 
analysis of H493 is that r[Nll] can be used for the determi- 
nation of ionic abundances of all the low-ionization species. 

From the electron temperatures in Table [3] the ionic abun- 
dances presented in Table |4] were derived. To calculate the 
total abundances of oxygen, nitrogen and sulfur, the follow- 
ing equations were used: 



0/H=(0 + + 0^)/H + 
N/0 = N + /0 + 



(2) 



(3) 
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TABLE 2 

Reddening-corrected line fluxes 



Line 




H493 


H507 


H972 


H1013 


3727 


rn tti 


84 ± 4 


84 it 4 


1 3Q _L 7 
1 J7 It / 


221 ±11 


J 1 JU 


i_| i 9 
ri iz 




1 4-09 
1.0 It U.Z 




a n 4- 9 

J.U It U.Z 


■3,771 

J / / 1 


i_i i i 
ri 1 1 




A 7 4- 3 
<+. / It U.J 




J.J It U.Z 


"X1Q9. 
J / Vo 


i_i i n 

ri iu 


ft a _i_ a 

D.J It U.J 


S ^ 4- ^ 
J.J It U.J 




4.4 It U.Z 


'3.899 
joZZ 


rlc 1 








79 4- 0£ 
U. /Z It U.UO 


JOJJ 


n y 


O.J It U.4 


7 1 4-0 4 
/ . 1 It U.4- 




7 3, 4_ % 

1.3 It U.J 


3.8^8 
joOO 


LrNc 111 J 








a 1 4-0 9 
J. 1 it U.Z 


3,QftQ 


I_J o 

n o 


1 ft a _i_ 7 

lO.UxU. / 


1 9 < _|_ £ 
1Z.J it U.O 


1Q1 4-0 7 
lo. 1 it U. / 


1 £ 4 4- 7 
10.4 it U. / 


4026 


He I 








1 S4 4- 08 
1 . J4 It U.UO 


4072 t 


r<N tti 








1 27 ± 07 


A1 01 
4 1U1 


no 


9ft A 4- 1 1 

zo.'t in 1.1 


9£ A 4- 1 1 
ZO.t It 1 . 1 


97 ^ 4- 1 1 
Z / . J It 1.1 


lf\ 9 4-11 
ZO.Z It 1.1 


A 1 /l /I 
4 144 


rlc 1 








A 1 1 4- O Ol 
U. 1 J It U.UJ 


A9£7 
4ZO / 


P TT 
V. 11 








9S 4- (YX 
U.ZJ It U.UJ 


4 J4U 


H7 


Aft O. _|_ 1 C 
4U.J IT l.O 


46 4 it 1 9 


46 4 ± 1 7 


Af. 4 4-10 
4U.4 It 1.7 


4 303 \ 


LU inj 








94 4- (Y\ 
U.Z4 It U.UJ 


418R 

4JOO 


He I 








AO 4- OS 
U.4U It U.UJ 


AA71 
44 / 1 


LI., T 

rlc 1 


19-1-09 
1 .Z It U.Z 


9 9 4-09 
Z.Z It U.Z 


19 4-01 
1 .Z It U. 1 


^04-09 
3.y It U.Z 


46SR 


[Fe III] 








40 4- OS 
U.4U It U.UJ 


471 1 
4/11 


[Ar I v] 








3C 4_ 04 
U.JJ It U.U4 


AR£1 


rr a 
tip 


1 on -i- a 

1UU It 4 


1 00 4- A 
1UU It 4 


1 00 4- A 
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7236 


Cn 




1.02 ±0.18 




0.27 ± 0.03 


7281 


He I 








0.51 ±0.05 


7325 f 


[Oil] 




0.9 ± 0.2 




2.4 ± 0.2 


7751 


[Arm] 








1.54 ±0.11 


8545 


P15 








0.69 ± 0.06 


8598 


P14 








0.66 ± 0.07 


8665 


P13 








0.82 ± 0.09 


8750 


P12 








1.07 ±0.11 


8863 


Pll 




1.14 ±0.13 




1.35 ±0.12 


9015 


P10 


1.51 ±0.24 


1.39 ±0.15 




1.86 ±0.15 


9069 


[Sin] 


10 ±2 


9±2 


17 ±3 


24 ±5 


C(W) 




0.48 ± 0.05 


0.74 ± 0.05 


0.40 ± 0.05 


0.35 ± 0.05 


EW(H/3) 




59 ± 1 


48 ± 1 


34 ± 1 


146 ± 1 


F(H/3) a 




29.0 ± 1 


63.4 ± 1 


31.8 ± 1 


157.0 ± 1 



NOTE. — Some lines belonging to a doublet/multiplet are identified by a sin- 
gle wavelength, for example [O II] AA3726, 3729 (here [On]A3727), [S II] AA4069, 4076 
([SII]A4072) and [O II] AA7319.7320 + AA7329, 7330 ([O II] A7325). Auroral lines are 
identified by the t symbol. 

"Extinction-corrected total flux, in units of 10~ 15 erg s~' cm -2 . 
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TABLE 3 

Measured electron temperatures (K) 



H493 



H1013 



1LWU,J 4959+5007 
T[OlI] ^ 

IT ] 6717+6731 
T [ Sln ] 9069 3 +9 2 532 



7700 ± 250 

7690 ± 200 
5960 ± 340 8360 ± 150 

8160 ±340 
7330 ± 590 7680 ± 160 



TABLE 4 

Ionic and total abundances from 
collision ally excited lines 





H493 


H1013 


0+/H+ 


5.4 ±2.3 x 10 -4 


2.3 ± 0.2 x 10"* 


0++/H+ 


1.4±0.8x 1CT 5 


9.9 ± 1.4 x 10~ 5 


N+/H+ 


1.0±0.2x 1CT 4 


2.6±0.1 x 10" 5 


S+/H+ 


6.1 ± 1.4 x 10" 6 


1.2±0.1 x 10~ 6 


S++/H+ 


3.1 ±0.7 x 10~ 6 


6.7 ± 0.3 x 10~ 6 


Ar++/H + 


3.3 ± 0.8 x 10~ 7 


1.3±0.1 x 10~ 6 


Ne++/H+ 




1.1 ±0.2 x 10~ 5 


12 + log(0/H). 


8.74 ±0.17 


8.52±0.05 


log(N/0) 


-0.72 ±0.22 


-0.95 ±0.05 


log(S/0) 


-1.78 ±0.20 


-1.61 ±0.05 



as suggested by Peimbert & Costero ( 1969), and 



S+ + S++ 



1-1 



O+v 

~o 



l/a 



(4) 



following Stasi riskal (119781) with a = 2.5. This last equa- 
tion yields very small corrections (~ 1% or less) to the sulfur 
abundance in both H493 and H1013, due to the low amount 
of S 3+ expected in l ow-excitation nebulae (see also KBG03, 
iBresolin et all 2004). The O/H, N/O and S/O abundance ra- 
tios thus derived are included in Table0] 

KBG03 determined the oxygen abundance gradient in 
M101 from 20 Hll regions in which the electron tem- 
perature of the gas could be determined from ratios of 
auroral lines to nebular lines. The galactocentric distances 
of these Hll regions vary between R/Rq = 0.19 (H1013) 
and R/Rq = 1 .25 (SDH323). The new measurements pre- 
sented here add one further point closer to the center of the 
galaxy, H493 at R/Rq = 0.10, and allow a determination of 
the O" 1 " 1 " abundance of H1013 from a direct measurement of 
[O III] A4363. The resulting O/H abundance is ~0.2 dex lower 
than measured by KBG03: 12 + log (O/H) = 8.52 ±0.06 
(this work) vs. 12 + log(0/H) = 8.71 ±0.05 (KBG03). 
The differences in the adopted + temperature (8070 K 
vs. 7600 K) and in the 0++ temperature (7700 K vs. 6600 K) 
are responsible, in approximately equal proportions, for 
the discrepancy in abundance. The O" 1 " 1 " zone temper- 
ature was derived by KBG03 from r[Sm], measured 
from the [S III] A6312/( A9069+ A9532) line ratio, and 
the application of Equation (1). The r[Sm] value de- 
termined in the current work is ^500 K higher than the 




FIG. 2. — The oxygen abundance gradient in M101 in terms of the fractional 
isophotal radius, R/Rq, determined from the KBG03 measurements (open 
circles) and the inner regions H493 and H1013 studied here (full dots). All 
the data points represent H II regions where the electron temperature has been 
measured directly from auroral-to-nebular line ratios. In all cases, except for 
the innermost region, H493, and H336 (= Searle 5, point at R/Rq = 0.22) the 
[O III] temperature has been derived from the A4363 line. In this plot the 
straight line represents the weighted linear least square fit, given by Eq. (5). 



KBG03 value, as a result of a ~35% smaller flux in 
the [S III] AA9069, 9532 lines. I point out that Equation 
(1), together with the value of T[S III] in Table [3] would 
provide T[0 III] =7200 ± 200 K, i.e. 500 K smaller than 
measured directly from JO ml A4363/( A4959 + A5007) . 
The relation proposed by Perez-Montero & Diaz (2005), 
r[S m] = 1.05 T[0 in] - 800 K, would yield instead 
T[0 III] = 8100 ± 150 K. This result is a warning that 
relations derived from photoionization models, such as 
Equation (1), can work well in a statistical sense, but can fail 
for individual H II regions when uncertainties in T e smaller 
than a few hundred degrees are required. 

A weighted linear least square fit to the KBG03 data points, 
supplemented by the abundances determined here for H409 
and H1013, provides the following solution for the radial oxy- 
gen abundance gradient in M101: 

12 + log(0/H) = 8.75 (±0.05) -0.90 (±0.07) R/R . (5) 

This relation is virtually the same as the one determined by 
KBG03. The innermost region studied in this paper, H493, 
allows us to extend the empirical abundance gradient in M101 
to 1.5 arcmin from the galactic center. The oxygen abundance 
of H493, 12 + log(0/H) = 8.74±0.17, agrees well with the 
galactic abundance gradient, as shown in Fig. [2] 

4. HELIUM ABUNDANCES 

There are several He I recombination lines which could, in 
principle, be used to obtain the helium abundance, especially 
in the case of H1013 (see Table 2). However, many of these 
lines are faint, with small equivalent widths, and the abun- 
dances derived from them would be quite sensitive to even 
small corrections for underlying absorption. Some of the 
stronger lines are not well suited for the measurement of the 
helium abundance either, because they are partially blended 
with strong metal lines, as in the case of HelA5016 with 
[O III] A5007, or they can be strongly affected by optical depth 
effects (He I A7065). Therefore, only the two lines He I A5876 
and A6678 have been used. 

The He + /H + ionic ratios have been obtained from 
H eiA5876/H/3 and HeiA 6678/H^, adopt ing the emissivities 
of Beniar rmTetai] {l999) for He I and of IStorev & Hum mer 
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TABLE 5 

Ionic and total abundances from 
recombination lines 



H493 H1013 



He+/H+ 0.048 ±0.003 0.083 ±0.001 

He/H 0.10±0.01 0.10±0.01 

C++/H+ •■■ 2.3 ±0.3x10"" 

12+log(C++/H + ) •■■ 8.36±0.06 

12+log(C/H) •■• 8.66±0.07 



i 19951) for H I, both calculated at the T[0 III] temperature de- 
termined in the previous section. The average ionic abun- 
dances thus derived from A5876 and A6678 are reported in 
Table |3 In the case of the low excitation Hll region H493 I 
find He + /H + = 0.048, almost half the value found in H1013, 
He + /H + = 0.083. These helium ionic fractions should be 
compared to oxygen ionic fractions + /(0 + ±0 ++ ) of 0.97 
(H493) and 0.70 (H1013; a value of 0.55 is however derived 
when accounting for the presence of temperature fluctuations, 
as shown in § 6). 

A large fraction of the helium present in low-excitation neb- 
ulae is in the neutral form. To estimate its contribution to the 
total helium abundance, a grid of p hotoionization mode ls has 
been obtained with cloudy 94.00 dFerland et alll998t) . The 
helium ionization correction factor, ICF(He), has been esti- 
mated from the predicted relation between He / (He°±He + ) 
and + /(0 + ±0 ++ ) (He ++ is negligible at these low excita- 
tion levels). I found ICF(He) = 2.13 ±0.23 for H493 and 
ICF(He) = 1.28 ±0.09 for H1013(1.20±0.09 accounting for 
temperature fluctuations). The resulting helium abundances 
forH493 (He/H = 0.10±0.01) and H1013 [He/H = 0.11 ± 
0.01, or He/H = 0.10±0.01 if using + /(0 + ±0 ++ ) =0.55] 
are in agreement, within the uncertainties. 

5. THE CARBON ABUNDANCE OF H1013 

The C II A4267 recombination line has been detected in 
H1013, the brightest of the H II regions studied in this work 
(Fig-EJ- Esteba net al.l J2002) measured this line, as well as 
various O II lines, in two H II regions of lower metallicity in 
M101, NGC 5471 and NGC 5461. Metal recombination lines 
allow the determination of ionic abundances that are virtu- 
ally independent of the electron temperature, thanks to the 
fact that they have a similar temperature dependence as the 
hydrogen recombination lines, and therefore allow very im- 
portant checks on the abundances obtained from collisionally 
excited lines. Unfortunately, no O II lines could be measured 
in the LRIS spectra. This is likely a result of the low spec- 
tral resolution, combined with the fact that the strongest oxy- 
gen recombination lines around 4650 A fall in the spectral 
region occupied by the strong W-R blue bump. Besides, at 
the low excitation levels typical of the H II regions considered 
here, most of the oxygen is in the + form rather than ++ , 
as shown in the previous section. As a consequence, the O II 
recombination lines will be very faint. This problem affects 
most metal-rich H II regions, which are generally of low exci- 
tation (see Sect. 7). 

Because of the lack of collisionally excited lines from car- 
bon in the optical spectra of Hll regions, and the current 
unavailability of UV spectrographs working in space for the 
study of lines like C III] A1909, the derivation of the C ++ /H + 



ratio and of the total C/H abundance in extragalactic nebulae 
from the detection of recombination lines is very important to 
study the chemical evolution of spiral galaxies. In the case 
of the center of M101, it is interesting to test whether the 
C II A4267 recombination line provides a carbon abundance 
that is consistent with the measured oxygen abundance. For 
this purpose, I have adopted the C II eff ective recombination 
coefficient c<cii from lDavev et al.l J2000) and used the follow- 
ing expression to derive the abundance of C ++ : 

= 1(4267) 4267 a m 
H+ I(H/3) 4861 acn 

where the H/3 effective recombination coefficient a hp has 
been obtained from the IStorev & Hummeil i ll 9951) emissivi- 
ties. Both a's have been calculated at the r[Olll] temper- 
ature. The resulting ionic abundance fraction is C ++ /H + = 
2.3 x 10" 4 , equivalent to 12±log(C ++ /H + ) = 8.36 ±0.06 (Ta- 
bleEJ. 

In order to account for the C + contribution to the total C 
abundance (C 3+ is negligible in the low-excitation regime), 
a procedure similar to the one adopted for the helium ICF 
has been followed, i.e. the C + /C fraction as a function of 
O" 1 "/ O predicted by photoionization models was considered. 
The result C + /C = 0.61 ± 0.02 has been obtained. The 
models by iGarnett et alJ ( 119991) would yield the same re- 
sult. Adopting instead + /0 = 0.55, as determined in § 6, 
I obtain instead C + /C = 0.50 ± 0.02. The latter value is 
adopted, to infer a total carbon abundance for H1013 of 
12+log(C/H) = 8.66 ±0.06. This value corresponds to 1.9 
times the solar value ( Asnlund et a 1. 2005). Fig. |4] compares 
the O/H and C/O abundance determined for H1013 with data 
from the literature. Abundances obtained from collisionally 
excit ed lines (Hll regions in spiral and irregular ga laxies 
from Garn ett et alJl^fl997l[T999tlKobulnickv & Skillman 
1998) and from recombin ation lines (H II reg i ons in the Milky 
Way and other spirals. lEsteban et al J 120021 l2005t |Peimbert 
l2003tlPeimbert et alJl2005t) are shown. The H1013 datapoint 
should be compared with the latter set only, but the collision- 
ally excited line results are also included for completeness 
(the dotted lines connect H II regions in common between the 
two sets: NGC 5461 and NGC 5471 in M101, and 30 Doradus 
in the LMC). 

6. TEMPERATURE FLUCTUATIONS IN H1013 

Thanks to the excellent blue sensitivity of LRIS, the Balmer 
discontinuity at 3646 A is very well observed in the nebular 
spectra presented in this work (the case of H1013 can be seen 
in the top panel of Fig. [3). The nebular Balmer jump offers 
the opportunity to derive an electron temperature that is inde- 
pendent of the collisionally excited lines of metals, and rep- 
resentative of the hydrogen continuum s pectrum, T(Bac). As 
shown originally by Peimbert ( 1967) and Peimber t & Costerol 
( 1969), the comparison between two independently derived 
temperatures, e.g. T(Bac) and r[Olll], provides an indica- 
tion about the presence of temperature fluctuations in Hll 
regions and planetary nebulae. In the following I analyze 
the Hll region H1013, where the contribution to the nebu- 
lar continuum by underlying stars (direct and scattered) ap- 
pears to be negligible, as judged by the small equivalent width 
(0. 1 A) of the Balmer lines in absorption inferred during the 
redde ning correction. I follow the technique developed by 
Peimbert & Costero ( 1969) and summarized in several recent 
papers, for example bv lPeimbert et alJ fe004l) . 
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FIG. 3. — Portions of the spectrum of H 1013 showing the Balmer discontinuity region and the high-order Balmer lines (top), and the C II A4267 recombination 
line, the [O III] A4363 auroral line and the Wolf-Rayet broad emission features around 4650A (bottom). 
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FIG. 4. — C/O ratio as a function of the oxygen abundance O/H from differ- 
ent samples of H II regions. Small triangles represent nebulae in spirals and 
irregulars for which the abundances have been obtained from collisionally 
excited lines i Garn ett et alll99alT991 1999; Kobulnickv & Skillm arll99a . 
Recombination line abundances exist instead for the objects represented by 
open (extragalactic) and solid (Galactic) circles i Esteban et al. 2002, 2005; 
Peimbert 2003; Peimbert et al. 2005). The dotted lines connect H II regions 
in common between the two sets. The H1013 H II region is shown as a large 
open square. 



As a first step, the electron temperature from the Balmer 
jump has been derived. For this purpose, one needs to repro- 
duce the observed Balmer discontinuity, 7(3646") -7(3646 + ), 
accounting for the continuum processes at work: free-free 
and free-bound transitions of Hi, He I and He II, and two- 
photon decay of the 2 2 Sj/2 level of hydrogen. The corre- 
sponding temperature-dependent coefficients have been taken 
from lBrown & Mathews! i 19701) . neglecting the Hell contri- 
bution, because of the small excitation of H1013, and adopt- 
ing the He + /H + ratio in Table|5]for the calculation of the He I 



contribution. The spectral resolution of the HI 01 3 spectrum 
is not sufficient to allow a measurement of the real continuum 
in proximity of the Balmer discontinuity, because the blend- 
ing of the high-order Balmer lines creates a raised pseudo- 
continuum. I have therefore extrapolated a fit to the nebular 
continuum in the ~ 3800-3950 A wavelength range to the 
wavelength of the Balmer discontinuity. Using then the ratio 
of the discontinuity to the H/3 and H7 emission lines I have 
derived T (Bac) = 5000 ± 800 K, where the error is dominated 
by the extrapolation of the continuum. 

In order to derive the average temperature Tq and the mean 
square temperature fluctuation t 2 the following equations 
were used: 



r[oral=7b 



1+ (^-3)f 



(7) 



T(Bac) = 7b(l-1.67f 2 ) 



(8) 

yielding T = 5500 ± 700 K, t 2 = 0.06 ± 0.02. Comparable 
values of f 2 have been found for other extragalactic H II re- 
gions. For example, f 2 = 0.041 in NGC 5461 and t 2 = 0.074 
in NG C 5471 (both H II r egions located in M101 and stud- 
ied bvlEsteban et al 120021), t 2 = 0.064-0.098 in NGC 2363 
JGonz alez-Dekado et al. 1994 k t 2 = 0.0 54-0.076 in regions 
X and V of NGC 6822 jPeimbert et alJl2005l) . f 2 = 0.033 in 
30 Doradus (Peimb eri2003l) . 

The effect of temperature fluctuations is to increase the 
chemical abundance measured from collisionally excited lines 
determined assuming t 2 = 0. The ionic abundances in Ta- 
ble|4|were ther efore correc ted for the case t 2 ^ 0, following 
Peimbert & Costero ( 1969). The results are summarized in 
Table [6] As the table shows, the correction to the oxygen 
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TABLE 6 

H1013: IONIC AND TOTAL ABUNDANCES CORRECTED FOR 



t 2 = 0.00 t 2 = 0.06 



0+/H+ 2.3 ±0.2x10^ 4.7 ±0.4x10^ 

0++/H+ 9.9±1.4xl0~ 5 3.9 ±0.6x10^ 

N+/H+ 2.6-0.1 ID ' 5.1 • 0.2 :o 5 

S+/H+ 1.2±0.1xl0- 6 2.4±0.2xl0~ 6 

S++/H+ 6.7±0.3xl0^ 2.9±0.1xl0~ 5 

Ar++/H+ 1.3±0.1xl0 -6 4.1 ±0.3x10^ 

Ne++/H+ I.l±0.2xl0~ 5 4.7±0.8xl0~ 5 

12 + log(0/H) 8.52 ±0.05 8.93 ±0.05 

log(N/0) -0.95 ±0.05 -0.96 ±0.05 

log(S/0) -1.61 ±0.05 -1.47±0.07 



abundance amounts to about 0.4 dex. This upward correc- 
tions to O/H is somewhat larger than typically found compar- 
ing results from collisionally excited lines and recombination 
lines. The average correction in a numbe r of H II regions stud- 
ied in the Milky Way is 0.2 dex (see lGar cfa-Roias et al.l2006l 
Peimb ert & Peim bert 2005, and references therein), ranging 
between 0.11 dex and 0.25 dex. Slightly larger values are 
found in the extrag alactic Hll regions NGC 5461 (0.25 dex, 
Esteb an et alJIgOof) and Hubble V in NGC 6822 (0.29 dex, 
IPeimbert et al.l2005l) . 

In conclusion, accounting for the presence of tempera- 
ture fluctuations, the oxygen abundance of H1013 is 12 + 
log(0/H) = 8.93 ±0.06. Using the carbon abundance from 
§ 5, the carbon-to-oxygen ratio is log(C/0) = -0.27 ±0.11. 
The measurement of the Oil recombination lines in H1013 
would improve the accuracy of these determinations, and will 
require spectroscopic observations at higher spectral resolu- 
tion than the one used for the current work. 

7. DISCUSSION 

In this section I consider the results obtained for the inner 
Hll regions of M101 within the context of the calibration of 
nebular metallicity indicators based on strong lines, in par- 
ticular i?23- In general most of the following considerations 
could be applied to other strong-line diagnostics, some of 
which do not suffer some of the problems that affect /?23- The 
major one is the well-known double valuedness of the param- 
eter, that is the fact that a given R23 value corresponds to two 
values of the oxygen abundance. An emission line diagnos- 
tic that is monotonically increasing with abundance, such as 
[N ll]/[0 II], can be used to remove the degeneracy and place 
an H II region in the correct branch. For the metal-rich H II 
regions that are of interest here, only the upper branch cal- 
ibration is rele vant. S econdly, is sensiti ve to the ioniza- 
tion parameter ( Pere z-Montero & DfaJ2005T). This can be ac- 
counted for, in principle, using \Q ni]/\ Q III (M cGaug h[l99Tl) 
or [O lll]/([0 III] + [O II]) (the P index, lPilvuginll2001l) . De 
spite these drawbacks, R23 is preferentially considered in what 
follows because of its wide application in the literature, and 
because much of the recent empirical work deals with this in- 
dicator. 

7.1. Abundance biases 

This work has extended the auroral line method of abun- 
dance determination to a fractional galactocentric distance in 



M101 of R/Rq = 0.10 (~ 3 kpc) and to an oxygen abundance 
12 + log(0/H) = 8.74. The innermost Hll region studied by 
KBG03 was H1013 (R/R = 0.19). The oxygen abundance 
they determined from the [S III] and [N II] auroral lines has 
been revised to a ~0.2 dex lower value with the detection of 
[Om] A4363 in this paper. Both H1013 and the Hll region 
closest to the center analyzed here, H493, have been found to 
fit the oxygen abundance gradient in M101 quite well. A sys- 
tematic offset from this gradient would be expected if metal- 
rich H II regions are affected by strong abundance biases. The 
result obtained for H493 and H1013 then suggests that the 
biases due to large-scale temperature fluctuations in H II re- 
gions of approximately solar metallicity, 12 + log(0/H) = 8.7, 
are not as large as the maximum value predicted by photoion- 
ization models (a ^0.2 dex systematic effect in the case of 
abundances based on [O III] A4363, Stasinskal l2005l) . or that 
the expected abundance biases become relevant at a some- 
what larger metallicity, perhaps around 12 + log(O/H) = 9.0. 
In either case, the abundances obtained from au roral l ines in 
this and in previous papers (e.g. Bresolin et al. 2004) up to 
approximately the solar metallicity, i.e. what one finds in the 
very central regions of metal-rich spiral galaxies, appe ar to be 
robus t (a similar argument is also presented by Pilvugin et al. 
120061) . 

7.2. Metallicity dependent excitation 

The direct abundances measure d in recent year s from 
auroral-to -nebular line ratio s JCaste llanos et all 120021 
KBG03, iBresolin et all l200l 120051) have provided the 
necessary input for the empirical calibration of strong-line 
abundance i ndicators for metal-r ich H II regions. The 
P-method of PilvugirJ (120001 12001[) has been recently recali- 
brated at high metallicity by Pilvugin & Thuan (2005) on the 
basis of the recent r e -based abundances. As mentioned ear- 
lier, this empirical calibration leads to H II region abundances 
that are, in the case of objects located in the metal-rich, 
upper branch of R23, factors of 2-3 lower than those derived 
from calibrations based on photoionization models. In 
order to calibrate indicators that can provide abundances of 
star-forming galaxies across the universe with accuracies 
of, say, 0.1 dex, it is crucial to understand the origin of 
the discrepancies between different calibrations, and try to 
test them against additional empirical methods. The main 
difficulty is to provide accurate metallicities for H II regions 
in the log7?23 < 0.5 regime, corresponding approximately to 
(O/H) > (O/H) , due to the weakness of the direct abundance 
diagnostics. 

In Pilyugin's method, R23 depends not only 
on the oxygen abundance, but also on a mea- 
sure of the nebular ionization, parameterized by 
the excitation parameter P = [O III] AA4959, 5007 / 

(IP ml AA4959, 5007 + \ Q ill A3727) . This follows the 

original suggestion by McGaughf ( 119911) to introduce a 
two-dimensional characterization of nebular spectra. It 
is worth pointing out here the empirical finding that the 
extragalactic H II regions that extend along the upper branch 
of the R23 vs. (O/H) diagram below log/?23 = 0.5 tend to 
have a low excitation parameter. In fact, an increase in 
the metallicity (decreasing R23) corresponds, in a statistical 
sense, to a decrease in P, i.e. high metallicity H II regions 
with a measured value of T e tend to be also low excitation 
nebulae. This can be seen by plotting R23 against O/H, 
togeth er with curves of con stant P (see, for example, Fig. 12 
of lPilvugin & Thuanl2005l) . and can be viewed in light of the 
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FIG. 5. — Trend of decreasing equivalent width of H/3 with decreasing 
excitation parameter P. The d ots represent Hn regions samples in M101 
(KBG03, crosses), M51 i Bresolin et al. 2004, open circles), and in five addi- 
tional spiral galaxies (Bresolin et al. 2005, open triangles). The four M101 
H II regions studied in this work are labeled (square-dot symbols). 



correlation between galactocentric distance and equivalent 
width of the H/3 nebular emission line, EW( H/J), noted in 
the early abundance stu dies of spiral galaxies (Searle1 IT97lL 
LShields & Tinslevl 119761) . This is shown in Fig. |5] where 
log EW(H/ 3) is plotted against P for th e H II region samples 
of KBG03, Bre solin et all 120041 12005) and for the M101 H II 
regions studied in this paper. This plot illustrates the fact that 
the ionizing clusters of H II regions located in the metal-rich 
inner regions of spiral galaxies are characterized by smaller 
effective tempe ratures than those further away from the 
galactic nuclei ( Bresolin & Kennicutt 2002). A metallicity 
dependence of the ionization parameter U = Q H o /4-nR^nc 
(the ratio of ionizing photon density to atom density; R$ is 
the Stromgren radius), as well as the relatively narrow range 
in the observed U values, have been noted in stud ies of extra- 
galactic H II regions several times (e.g. [D opita & Evans 1986; 
Bresolin et al. 1999; Kewlev & Dopita 2002), and have been 



confirmed by recent investigations of the integ r ated s pectra 
of star-forming galaxi es JNagao et al.l I2Q06; Mai er et aD 
120061 ILlang et al.ll2006l) . The lower excitation of metal-rich 
nebulae can be explained by a combination of factors. 
First of all, stellar atmospheres of massive O stars become 
cooler with increasing me tallicity, as a result of enhanced 
line and win d blanketing (Masse v" et alJ 1200*51) . Moreover, 
Dopit aeTaTl (|2£>06b) have shown that the ionization param- 
eter is determined by the ratio of the ionizing photons (a 
decreasing function of metallicity, due to enhanced stellar 
atmosphere opacities) to the mechanical energy input from 
the ionizing O stars (an increasing function of metallicity, 
since the photon momentum is more efficiently transferred 
to the stellar wind as the ab undance of absorbers increases, 
iBresolin & Kudritzkil 12004). The resulting metallicity 
dependence derived by Dopita et al., given approximately 
by U oc Z~° 8 , is in good agree ment with some empirical 
measurements. For example, Bresolin et al. ( 1999) estimated 
a factor 4 decrease in U varying the metallicity from 0.2 Z e 
to l.OZ . 



7.3. Temperature fluctuations 

The analysis based on the classic auroral line method pro- 
vides lower limits to the oxygen abundance in H II regions. 
Accounting for the presence of temperature fluctuations leads 
to higher abundances relative to what is derived by assum- 
ing a constant temperature. Temperature fluctuations in neb- 
ulae can be estimated in different ways. As done here in 
the case of H1013, one can compare the temperature derived 
from the Balmer (or Paschen) continuum with r[Olll] (or 
other ionic temperature measured from an auroral-to-nebular 
line ratio). In virtually all of the Galactic and extragalactic 
Hll regions with 12 + (0/H)>8.1 in which this comparison 
has been carried out, the for mer is found to be si g nificantly 
lower than the latter (e.g. [P eimbert et al. 2000; Peimbert 
l2003UGarcfa-Roias et al.ll2006l) . For a sample of metal-poor 
[12 + log(0/H) = 7.1-8.3] emission-line galaxies lGuseva et ail 
(2006) found, instea d, that the two temp eratures do not differ 
significantly, while Hagele et al. (2006) found that only one 
out of three metal-poor [12 + log(0/H) = 7. 9-8.0] Hll galax- 
ies shows significant temperature fluctuations. 

It is also possible to estimate t 2 by requiring that the abun- 
dances obtained from col lisionally e xcited lines and from re- 
combination lines match (Peimbert et al. 199 31), and from the 
analy sis of helium recombination lines (Peimbert et alJl2000l 
2005). These methods have provided consistent results when- 
ever a comparison has b een attempted ( Esteban et alJl2004l 
Garcfa-Roias et al. 2004). In extragalactic H II regions of high 
metal content the most feasible approach is the one adopted 
here for H1013, since metal recombination lines are difficult 
to detect. As the excitation of the H II regions decreases with 
increasing metallicity (Fig. the O II lines become fainter, 
despite the increase in the oxygen abundance. Even so, further 
efforts to detect metal recombination lines in high-metallicity 
nebulae, similar to what has been done for C II A4267 in 
H1013, would be very valuable. A somewhat higher spec- 
tral resolution (R ~ 5,000- 10,000) than the one used in this 
work on M101 will aid in the detection of these lines, as well 
as allow to disentangle lines belonging to a multiplet, such as 
Oil. 

The effects of temperature fluctuations on the derivation of 
nebular abundances need to be accounted for in the calibration 
of metallicity indicators based on strong lines. In Galactic and 
extragalactic H II regions an increase in oxygen abundance 
by a typical factor of 2 to 3 is found relative to the results 
obtained from collisionall y excited lines under t he assump- 
tion of t 2 = 0. Recently, Peimb ert et alJ J2P06) have com- 
piled a list of H II regions where the abundances have been 
measured, mostly by the same authors, via metal recombi- 
nation lines, as well as from collisionally excited lines, with 
the goal of establishing a preliminary calibration of R23 as 
a function of (O/H) from the recombination line technique. 
Fig- El shows the abundances measured in H1013 from this 
paper, for t 2 — (large open square) and for t 2 = 0.06 (large 
full square), i n the 12 + log(Q / H) vs . log R23 plane. The H II 
regions in the Peim bert et alJ J2006I) compilation are shown 
by the small open squares (abundances from collisionally ex- 
cited lines, t 2 = 0) and by the small filled squares (abun- 
dances from metal recombinations lines or from an estimate 
of t 2 based on the Balmer continuum temperature). All ob- 
jects plotted in the diagram belong to the upper branch of R23, 
including those in the compa rison sample of metal- rich H II 
regions from KBG03 and Bresolin et al. ( 2004 120051 crosses, 
open triangles and open circles, respectively). The H II region 
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H1013inM101 is currently the extragalactic nebula that lies 
at the lowest R23 value among those that currently have a t 2 
estimate. Its oxygen and carbon abundances are among the 
highest measured so far in this relatively small sample (18 
objects). 

In Fig. |6] the empirical measurements are compar ed with 
the P-method calibration of Pilvugi n~& Thuanl (120 05) for two 
representative values of P (0.2 and 0.6, solid curves). The 
nebular line-based abundances, derived by neglecting tem- 
perature fluctuations (open symbols), are well represented by 
these curves, and show the tendency mentioned above for a 
decreasing excitation parameter as the metallicity increases. 
The abundances derived from metal-recombination lines or 
accounting for temperature inhomogeneities (solid squares) 
are displaced upward by about 0.3 dex, approaching the theo- 
retical calculati ons from phot oionization models, taken from 
iKobulnickv & Kewlevl (|2p04) for two different values of the 
ionization parameter, 10 7 cms _1 and 10 8 cms _1 (these cor- 
respond to values of the dimensionless ionization parameter 
U = q/c of approximately log U = -3.5 and log U = -2.5, 
respectively). Despite the difficulties that photoionization 
models have at reproducing emission lines of real nebulae 
that are very sensitive to T e , such as [Olll]A4363, they are 
more successful with lines that have a smaller T e dependence, 
such as [O II] A3727 and [O III] A5007. This would explain 
the relatively good agreement found in Fig. 5 with the abun- 
dance determin ations that account for temperature fluctua- 
tions |Peimbert et al. 2006). Still, the problem remains that 
standard photoionization models predict values for t 2 in H II 
regions that are very small, between 0.001 and 0. 01, com- 
pared to observed values between 0.02 and 0.09 (Peimbert 
ll995llPerezlT997i) . 

In future observations it will be important to extend the 
sample of objects where the abundances are obtained from re- 
combination lines to lower values of R23, down to log R23 ~ 0, 
if we want to establish a calibration of R23 based on recom- 
bination lines and/or t 2 determinations that can be used at 
the highest metallicities. This is made somewhat problem- 
atic by the weakness of the metal recombination lines and 
by the difficulty in obtaining reliable B aimer continuum tem- 
peratures. Therefore, it is still valuable to measure nebular 
abundances from [O III] A4363 and other auroral lines in the 
high-metallicity regime, since, under the temperature fluctua- 
tion paradigm, they appear to provide good estimates for the 
lower limits in abundance. 

7.4. [Nn]/[On] 

As mentioned before, the [N ll]/[0 II] ratio can be used to 
remove the ambiguit y between t he low er and upper branches 
of /?23- In fact, Dopitaetal. (2000) showed that this ra- 
tio is an excellent abundance indicator by itself, superior to 
R23 because of its mono tonic dependence on metallicity and 
its insensitivity to the nebular excitation. In Fig. 01 show 
x = log([Nn] A6583/[0 11] A3727) as a function of 7>based 
oxygen abundances for the same extragalactic H II regions in- 
cluded in Fig. |6] supplement ed at the l ower -abundan ce end 
by the sam ples studied by|Garnett et al.l ill 997T). Ivan Zee et al. 
(1998) and lvan Zee & Havnesl ( 120061) dPerez-Montero & Diaz 
2005 showed a similar plot that includes a larger sample of 
objects). A simple fit to these data (dashed line) is: 

12 + log(0/H) = 8.66 + 0.36x-0.17x 2 (9) 
According to the models bv lDopita et al.l ( 1200 6a. their Fig. 7), 
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FIG. 6. — The upper, high-metallicity branch of /?23- Sa mples of metal- 
rich H II regions are taken from KBG03 (crosses) and Bresolin et afl 120041 
2005, open triangles and open circles, respectively). Small squares represent 
Galactic and extragalactic H II regions, compiled by Peimbe rt et alj 1200fl) . 
where abundances have been derived both from collisionally excited lines 
under the assumption t 2 = (open symbols) and from metal recombination 
lines or an estimate of t 2 (full symbols). The Hll regions H493 and H1013 
studied in this paper are represented by the lar ge square symbols. The P- 
method calibration of Pilvuain & Thuan 1 2005) is shown for P = 0.2 and 
P = 0.6 (full lines) . The Rn ca libration based on photoionization models 
by Kobulnickv & Kewlev 1 2004 ) is drawn for two values of the ionization 
parameter, q = 10 7 cm s (log U = -3.5) and q = 10 s cm s (log U = -2.5) 
(dashed lines). 



a given x corresponds with very good approximation to a 
unique value of the metallicity. The 12 + log(0/H) values 
taken from the Dopita et al. models are indicated in the bot- 
tom portion of Fig. [7\ This helps to illustrate the typical re- 
sult that abundance diagnostics calibrated via photoionization 
models provide oxygen abundances about 0.2 dex larger than 
those obtained from direct methods. From what has been 
shown in Section 7.3, temperature fluctuations might offer an 
explanation for this discrepancy. The empirical calibration of 
the [N ll]/[0 II] abundance diagnostic also shows that, despite 
the fact that this line ratio is not sensitive to the ionization pa- 
rameter, the accuracy of the oxygen abundances that can be 
derived from it is not better than what can be obtained from 
R23. The rms of the fit shown in Fig.Qis 0.20 dex, while lim- 
iting the fit to objects that are on the R23 upper branch (the 
extragalactic Hll regions shown in Fig.|6j gives rms ~ 0.15, 
the same one obtains by fitting these abundances as a function 
of log R23. This suggests that our poor knowledge of addi- 
tional parameters affecting the spectra of H II regions, such 
as the temperature and ionization structures, currently limits 
the accuracy in determining abundances from strong-line di- 
agnostics. 



8. CONCLUSIONS 

In this paper I have presented new deep spectra obtained 
with Keck/LRIS of four Hll regions located in the central, 
metal-rich (O/H above solar) zone of the spiral galaxy M101. 
The main results obtained are summarized as following: 

- electron temperatures have been measured from auroral-to- 
nebular line ratios for two Hll regions, H493 and H1013. 
The classic analysis based on collisionally excited lines pro- 
vides oxygen abundances of 12 + log(0/H) = 8.74±0.17 and 
12 + log(0/H) = 8.52±0.05 for these two objects, respectively. 
These measurements extend the direct determination of the 
radial abundance gradient in M101 to about 3 kpc from the 
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center (R/R = 0.10). 

- the C II A4267 recombination line has been measured in 
H1013. Accounting for the unseen C + /H + fraction with 
the aid of photoionization models, a total carbon abundance 
12 + log(C/H) = 8.66 ± 0.07 (approximately 1.9 times the so- 
lar value) is found. 

- for H1013 the comparison between the electron tempera- 
ture derived from the Balmer jump, T(Bac) = 5000 K, and 



the electron temperature derived from [Olll] A4363/A5007, 
T[0 III] = 7700 K, leads to a mean square temperature fluctu- 
ation t 2 = 0.06 ± 0.02, with a volume averaged temperature 
T = 5500 ±700 K. 

- correcting the abundances derived from the collisionally ex- 
cited lines for the effect of temperature inhomogeneities in- 
creases the oxygen abundance in H1013 by almost 0.4 dex, to 
12 + log(0/H) = 8.93 ± 0.05. Combining this result with the 
carbon abundance derived from C II A4267 gives log(C/0) = 
-0.27 ±0.11. 

- these results indicate that there is still no empirical verifica- 
tion of the strong abundance biases from temperature gradi- 
ents predicted theoretically for metal-rich H II regions. 

- temperature fluctuations can quantitatively explain the cur- 
rent abundance discrepancy of ^0.2-0.3 dex between empiri- 
cal indicators and theoretical models. 

Future observational work needs to target high-metallicity, 
low excitation H II regions in the local Universe with medium 
resolution spectroscopy, in search for metal recombination 
lines. From these, chemical abundances that are virtually 
independent of the temperature structure of H II regions can 
be derived. This would allow more accurate calibrations of 
strong-line metallicity indicators. 



It is a pleasure to thank Manuel Peimbert for encourage- 
ment and suggestions. 
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